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i—i . Abstract 

■ 

■ Form factors parameterizing radiative leptonic decays of heavy 

. mesons (B + — > ^yl + ui) for photon energy are computed in the language 

Qh| of dispersion relation. The contributing states to the absorptive part 

pi ■ in the dispersion relation are the multiparticle continum, estimated by 

quark triangle graph and resonances with quantum numbers \~ and 
1 + which includes B* and B* A and thier radial excitations, which model 
' the higher state contributions. Constraints provided by the asymp- 

^ ■ totic behavior of the structure dependent amplitude, Ward Identities 

and gauge invariance are used to provide useful information for pa- 
rameters needed. The couplings qbb*^ and /ss^ 7 are predicted if we 
restrict to first radial excitation; otherwise using these as an input the 
radiative decay coupling constants for radial excitations are predicted. 
The value of the branching ratio for the process B + —* 7/i + z^ is found 
to be in the range 0.5 x 10 -6 . A detailed comparison is given with 
other approaches. 



1 Introduction 

In spite of small branching ratio, the radiative I?-meson decay (B — > Iw-f) 
is of viable interest because it contains important information about weak 
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and hadronic interactions of 5-meson. Furthermore, with the introduction of 
5-factories LHCb, BaBar, Belle and CLEOb, the radiative I?-meson decay 
can be studied with enough statistics. Preliminary data from the CLEO 
collaboration indicates the limit on the branching ratio B(B — > lwy) which 
is: 

B(B -> ez/ e 7) < 2.0 x 10" 4 
B(B -> nvfl) < 5.2 x 10~ 5 

at 90% confidence level [1]. With the better statistics expected from the 
upcoming B factories, the observation and experimental study of this decay 
could become soon feasible. It is therefore of some interest to have a good 
theoretical control over the theoretical uncertainties affecting the relevant 
matrix elements. 

The radiative leptonic decay B + — > l + v\ r ) has received a great deal of 
attention in the literature [2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13] as a means of 
probing aspects of the strong and weak interactions of a heavy quark system. 
The presence of the additional photon in the final state can compensate for 
the helicity suppression of the decay rate present in purely leptonic mode. As 
a result, the branching ratio for the radiative leptonic mode can be as large 
as 10~ 6 for the ji + case [10], which would open up a possibility for directly 
measuring the decay constant f B [7]. A study of this decay can offer also 
useful information about the Cabibbo-Kobayashi-Maskawa (CKM) matrix 
element \V u b\ [14, 15]. 

In the radiative 5-decay process, there are two contributions to the am- 
plitude: 

1. inner bremsstrahlung (IB) and 

2. the structure dependent (SD) contribution which depends on the vector 
and axial vector form factor Fy and Fa respectively. 

The IB contribution to the decay amplitude is associated with the tree 
diagrams shown in Figs, la and lb, and SD contribution is associated with 
Fig. lc. 

In this paper, we will study the radiative leptonic B decays of B + — > 
Z + z//7. The IB part is still helicity suppressed [2], while the SD one is free 
of the suppression [16]. Therefore, the radiative decay rates of B + — > / + z^7 
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(Z = e, //) could have an enhancement with respect to the purely leptonic 
modes of B + — > l + ui due to the SD contributions in spite of the electro- 
magnetic coupling constant a. With the possible large branching ratios, the 
radiative leptonic B decays could be measured in the future experiments at 
hadronic colliders, such as BTeV and CERN Large Hadron Collider (LHC-B) 
experiments [17]. 

The paper is organized as follows. In Sec. 2, we present the decay kine- 
matics and current matrix elements for B + — > l + u^. In section 3, we dis- 
cuss the various contributions to the absorptive part of the SD amplitude 
iH^ V) needed in the dispersion relation. This include multiparticle contin- 
uum and resonances with quantum numbers 1~ and 1 + . The resonances 
include B* and B* A mesons and their radial excitations, which model the 
higher states. The continuum is estimated by quark triangle graphs. In Sec. 
4, the asymptotic behavior of the SD amplitude is studied. This provides a 
usual constraint on the residues of the resonance contribution, in terms of 
the continuum contribution. In Sec. 5, we discuss Ward Identities which 
together with gauge invariance relates various form factors. These identities 
which are expected to hold below the resonance regime, fix the normalization 
of the forms at q 2 = in terms of universal function g + (0) as well as another 
constraint on the residues. Thus in our approach, a parametrization of q 2 
dependence of form factors is not approximated by single pole contribution. 
But this parametrization is dictated by considerations mentioned above and 
also predict the coupling constants of 1~ and 1 + resonances with photon if 
we restrict to one radial excitation; otherwise using these as input, the radia- 
tive coupling constants of radial excitations are predicted. In this and other 
aspects our approach is different from the others mentioned previously. Our 
approach is closest to the one used in [18] for B — > nlui. We calculate the 
decay branching ratios in Sec. 5. We give our conclusions in Sec. 6. 

2 Decay kinematics and current matrix ele- 
ments 

We consider the decay 

B +(p) - (1) 
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where / stands for e or /x, and 7 is a real photon with k 2 = 0. The decay 
amplitude for radiative leptonic decay of B + — > / + f/7 can be written in two 
parts, M IB and M SD , as follows: 

M (B + -> Z+z/,7) = M /B + M SD (2) 

in terms of two emission types of real photon from B + — > l + u t . They are 
given by [19, 20, 21, 22] 

M IB = te^=V ub f Bmi e;L» (3) 

M SD = -i^V ub f Bmi e;Hn u (4) 

with 

Z" = «(p,)7 M (l + 7 5 )^,^), (6) 

^ = iF v {q 2 )e^k aP p-F A {q 2 ){ p .kg^ - jTk v ) , (7) 

g M = ( p -kr = ( Pl +p u r. (8) 

Here e* denotes the polarization vector of the photon with /c M e* (k) = 0. p, 
pi, p u , and k are the four momenta of B + , Z + , z/, and 7, respectively, is the 
polarization vector of the l + , f B is the B meson decay constant, and F A) F v 
stand for two Lorentz invariant amplitudes (form factors). 

The term proportional to L M in (5) does not contain unknown quantities- 
it is determined by the amplitude of the non-radiative decay B + — > 
This part of the amplitude is usually referred as "inner bremsstrahlung (IB) 
contribution", whereas the term proportional to H^ u is called "structure 
dependent (SD) contribution". 

The form factor Fa (Fy) is related to the matrix element of the axial 
(vector) current. The factors f B and F V>A are defined by 

(0|«7" 7 5&|fl(p)) = -if B p^ (9) 
(7(A0|«7^|5(P)> = -{(e*-p)k»-^(p-k)}F A (q 2 ) (10) 
^{k^u^Bip)) = - ie ^elp a kpF v (q 2 ) (11) 
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In our phase convention, the form factors Fa and Fy are real in the 
physical region 

mf < q 2 < M 2 B 

where q is the momentum transfer. The kinematics of the decay needs two 
variables, for which we choose the conventional quantities and in the rest 
frame of B 

2p-k 2£ 7 

x = (12) 

lp-Pi 2E l 
V = ~Ml = W B (13) 

and the angle 6^ 7 between the photon and the charged lepton is related to x 
and y by 

1 (2 - y + V?/ 2 " 4n) (2-y- Vv 2 - 4r,) 

2; — i , ,. — — (14) 

2 2-y + Vj/ a -4r,cos0, 7 

In terms of these quantities, one can write the momentum transfer as 

q 2 = M 2 B (l-x), (k 2 = o). (15) 
We write the physical region of x and y as 

< x<l-ri, (16) 

+ < y < 1 + n, (17) 

1 — x 

where 

_ f 9.329 x 10- 9 (Z = e) , 
' ~ M 2 B ~ \ 4.005 x 10- 4 (l = fi). [ } 

3 Dispersion Relations 

The structure dependent part, H**" is given by 

= 1 J d 4 xe ik -* (0 \T UL(x)JZ(0))\B(p)) (19) 

We note that [23] 

ik^ = if BPv , (20) 
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so that for the real photon we can write 

H iw = jji»> + $3——— (21) 
p ■ k 

where k fl H flu = and H^ v is parametrized as in Eq.(7). The second term in 
(21) is absorbed in Mjb- The absorptive part is 

Abs[iH^\ = yd 4 xe ik -*(0\[j» m (x),J»(0)}\B(p)) 

= ^(2^) 4 [E (0 n) (n\r 2 (0)\ B(p)} 5\k - p n ) 

-E (0 |^(0)| n) (n |^(0)| B{p)) 5\k + p n - p)] (22) 

n 

The 5-function in the first term implies p 2 n = k 2 = and since there is 
no real particle with zero mass, the first term does not contribute. Thus 
contributing to the absorptive part are all possible intermediate states that 
couple to .87 and annihilated by the weak vertex (0 \J% (0)| n). These include 
the multiparticle continuum as well as resonances with quantum numbers 1~ 
and 1+ Thus (t = q 2 ) 



, fB* A B~r , 

FAt) ~ Mj7^t fB - + 



(23) 



The ellipses stand for contributions from higher states with the same quan- 
tum numbers. The couplings Qbb*^ and /b*b 7 are defined as 



B*-(q,r ] ) 1 (k,e)\B-(P)) = tg B *B J e ap ^e* a q p r ] *^ 

B* A -(q, V ) 1 (k,e)\B-(P)) = ig B ^{e*.rf)-if B ^{q.e*){k.rf) 
(0|m7"6|S*(g,J7)> = f B .rf 
(0\iuYl 5 b\B A (q,r))) = f B *y (24) 

We assume that the contributions from the radial excitations of B* and B* A 
dominate the higher state contribution. Thus we write 

f a) - Rv 1 v Rv < 1 1 r 2 gg!!M , s 1 

V[t) - l-t/M 2 B , 2 -^l-t/M 2 ,7cJs s-t-ie + "' 
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(25) 

FA{t) - l-t/Ml, + \l-t/Ml, + nJ So s-t- l e dS + --- 

A A i 

where ellipses stands for the contributions from the region for much larger 
than the physical mass of heavy resonances up to oo. Here, M is a cut off 
near the first radial excitation of M B * or M b * a and So = M B + m n , and 

„ 9bb*2 , 

Rv ~ ~Mi: fB * 

(26) 

Ry i and i?^ are the corresponding quantities for the radial excitations with 
masses M B * and M B * . In the next section we develop the constraints on 

i A^ 

some of the parameters appearing in the above equations. 

If we model the continuum contribution by quark triangular graph (sim- 
ilar calculations exist in the literature [24]), we obtain 

pCont pCont fa J Qu Qb n A x [ 1 / 07 x 

Fv ~ Fa -W b \-k-m- b {1 + mJ)i-?/mI (27) 



where 

together with the term 



A = M B — m b , (28) 



frt ^ x f P^P" f P^P V 
k ■ p k ■ p 

which appears in Eq. (21). As is well known (see for example Ref. [25]), the 
pole at q 2 = M B in Eq. (27) arises due to u (bj quark propagator which form 
one leg of quark A, the other legs are the part of B meson wave function. 



4 Asymptotic Behavior 

To get constraints on the residues Ri, it is useful to study the asymptotic 
behavior of form factors Fy and Fa- It has been argued that the behav- 
ior of form factor for very large values of \t\ can be estimated reliably in 
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perturbative QCD processes [pQCD][26, 27, 18]. For t and for \t\ much 
larger than the physical mass of heavy resonances, pQCD should yield a very 
good approximation to the form factors. First we note that by vector meson 
dominance 

( 7 (k, e* (k)) \vrf (1 - 75) b\ B (p)> ~ Q U A (p (k, e* (k)) (1 - 7s) 6| B (p)> , 

Trip 

(29) 

where f p , having dimensions of mass, is defined as 

(0|«7"u|p(fc,e(fc))) = ^-e' 4 (30) 

trip 

Then using the methods employed in [27], it is eaisy to calculate [only the 
diagram where gluon is emitted by the light quark in (bu) bound state and 
absorbed by the heavy quark contributes and is by itself gauge invariant] 

ppQCD. 



7 pQCD jppQCD 



Qu 



fp \ 32vra s (t) 



VTLpj 3 



(/*/„)"»* Q In e)l (31) 



Here 

e~o(^-) (32) 



m B J 

and is governed by the tail end of the B meson wave function characterized 
by e. 

Now the asymptotic behavior of Eq. (25), is given by 



(33) 



F(q 2 )^ RM 2 + Y R t M 2 + - / %F Cont (s)d 

v ' q 2 [ i 7T J So 

Since F P ® CD (t) is a reliable approximation to the form factor for t — > —00, 
and (tF pQCD ^j -> in this limit, it follows that 

RM 2 + £i2 i M? + c~0, (34) 



where we have defined 

r-M 2 



%F Cont (s)ds. (35) 

7C J So 



The convergence relation (34) is a model-independent result and constitutes 
a very binding constraint for model building. In other words, the various con- 
tributions in Eq. (33) may be individully much larger than thef tF pQCD (t)j 




due to a s (t) / 1 suppression, but there must be large cancellations among the 
non-perturbative contributions in (33). This is in the spirit of ref. [18]. We 
will explore the resonant contribution (in our model) in order to understand 
the effect of Eq. (34) on the behavior of form factors in the physical region. 
The imposition of this constraint will lead to a very distinct behavior of the 
photon momentum distribution, independently of how many resonances we 
choose to keep. As the radial excitations of B* become heavier, they are 
less relevant to the form factors since the spacing between the consecutive 
radial excitations are expected to become narrower and narrower [28]. Thus, 
heavier resonances contribute with a smaller value even in the narrow width 
approximation. Furthermore, as finite widths are considered, the contri- 
bution of heavier and thus broader excitations are additionally suppressed. 
This shows that the truncation of the sum over resonances is a reasonable 
approximation. 

For the resonances stated above we will study a constrained dispersive 
model where only the first two radial excitations are kept. This is mainly 
for the reason mentioned above. On the other hand, the "minimal" choice of 
keeping only one radial excitation will determine R\ in terms of R. The other 
necessary ingredient to specify the model is the knowledge of the spectrum 
of radial excitations. These resonances [(2S) and (3S) excitations of B*] have 
not yet been observed in the B systems. We will then rely on potential model 
calculations for their masses [28]. These models have been very successful in 
predicting the masses of orbitally excited states and as such we are confident 
that the position of the radial excitations does not introduce a sizeable un- 
certanity. The resultant spectrum explicitly shows that the spacing among 
IS, 2S, 3S states are, to leading order, independent of heavy quark mass and, 
therefore, constitutes the property of the light degrees of freedom. The spec- 
trum of radial excitations is given in Table 1, where the subindices 1 and 2 
correspond to the 2S and 3S excitation of the B*, etc. Thus the convergence 
condition (34) now reads 



This condition leaves two free parameters Ri and R2 in the model. This 
results in the correct scaling of form factors with the heavy meson mass. 




RM 2 + RiMf + R 2 M\ + c = 0, 



(36) 
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Solving Eq. (36) for R 2 and using in Eq. (33), we obtain 



RM 2 (Mf - M 2 ) i? x M 2 (Mf - M 2 ) 



l r M2 M 2 - s 



(M 2 - q 2 ) (Mf - q 2 ) (Mj 2 - q 2 ) (Mf - q 2 ) M| - q 2 n Js s - q 2 

(37) 

If we model the continuum contribution by quark triangle graph as given in 
Eq. (27), we obtain 



C^pCont 



if) 



RM 2 (Mf - M 2 ) iW 2 (Mf - M 2 



M 2 



M 2 



(M 2 - g 2 ) (Mf - q 2 ) (M 2 - g 2 ) (Mf - q 2 ) (Mf - q 2 ) (M 2 - g 2 ) 

(38) 

where in the heavy quark limit M B = M* B = M and 



c = f B M B 



% • CO , 

A VM B 



(39) 



5 Ward Identities Constraints 

It is useful to define 

( 1 (k,e)\ma^q u b\B(p)} = -ie^ ] e* v k a p pF^q 2 ) (40) 
( 1 (k,e)\ma^ l5 q u b\B(p)) = [(q ■ k) e* M — (e* • q) k^ 1 } F 3 (q 2 ) (41) 

Now we will make use of Ward Identities and gauge invariance principle to 
relate different form factors. 

Usually, the gauge invariance is implemented by means of the Ward 
Identities; another way, essentially the same, is to consider what happens 
if the polarization vector of an external (real) photon is replaced by its four- 
momentum. The result is zero, provided that one considers all diagrams 
where this particular photon is connected in all possible ways to a charge 
carrying line. In this way one understands the connection between gauge 
invariance and charge conservation. The Ward Identities 1 used to relate dif- 
ferent form factors appearing in our process are: 

( 1 (k,e)\uia^q u b\B(p)) = -(m b + m q ) ( 7 (k,e) \u^b\ B{p)) 

+ (^ + P) (7 (k, e) \ub\B(p)) 

1 See ref.[29] for a detailed derivation of these Ward Identities. 
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= -(m b + m q )( 1 (k,e)\u^b\B(p)) (42) 
(7 (k, e) \uia^ 5 q v b\ B(p)) = (m b - m q ) (7 (k, e) \u^ l5 b\ B(p)) 

+(p» + W) ( 7 (fc,e) |« 7 56|5(p)> 
= {m b -m q ){ 1 (k,e)\u^ l5 b\B(p)) (43) 

where the matrix elements (7 (k, e) \ub\ B{p)) and (7 (fc, e) |-u^x 5 6| B(p)) vanish 
for real photon due to gauge invariance. 

Using the Ward Identities in Eqs.(40) and (41), and comparing the coef- 
ficients, we obtain [p ■ k = q ■ k, e* ■ p = e* ■ q] 

Fv(q 2 ) = -^—F^q 2 ) (44) 
m b + m q 

F A (q 2 ) = —F 3 (q 2 ) (45) 

m b - m q 

The results given in Eqs.(44) and (45) are model independent because these 
are derived by using Ward Identities. 

In order to make use of Ward Identities to relate different form factors, 
we define 

(7 (k, e) \iua af3 b\ B(p)) = -ie a(3pa e* p (k) [(p + k) a g + + q a g_] 

-iq ■ e*(k)e al3pa (p + k) p q a h 
-i [q a £f3paTt* p (k)(p + k) a q T - a <-> 0\hi 
~i i(P + k) a e fipaT e* p (k) (p + k) a q T - a <-»• f3] h 2 . 

(46) 

Since we have a real photon, gauge invariance requires that if we replace 
e p {k) by k^, the matrix element should vanish. This requires 

g + + g- + 2(q-k)h = (47) 

From Dirac algebra 

(7^75 = -le^Vaf,, (48) 

we can write 

{j(k,e)\iwT^ 5 b\B{p)) 
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-2 £ 



fivafl 



(7(fc,e) \iucT a pb\ B(p)) 



(e*»k» - e*"*") [<?+ - <?_ - (Ml + g 2 ) h, - (3M 2 B - q 2 ) h 2 

+ (e'V - e*V) [g + + 9- + (M 2 B - q 2 ) (h + h 2 ) 
-2q ■ e* (h - tn - h 2 ) (p^k u - p u k^) 



(49) 



The gauge invariance, namely, replacing e M by k^, the matrix element should 
be zero, does not give any new relation other than (47). Using this relation 
and 2k ■ q = M B — q 2 , we get 

( 1 {k,e)\iwr' Jv lf ,b\B(p)) 
= (e^fc" - e*"^) [2#+ + (M| - g 2 ) (h-h x - h 2 ) - 2q 2 h x - 2M 2 B h 2 
-[2k-q (e^p" - e* V) + 2g • e* (p"fc" - p v k»)\ (h - h x - h 2 ) (50) 

Contrary to what is stated in some literature, the gauge invariance does allow 
a second tensor structure in addition to {&k v — e v k^). 
This gives 

( 1 {k,e)\iqa^q vlb b\B{p)) = 2 (g + - q 2 h - (M 2 B - q 2 ) h 2 ) 

x (q ■ ke*" (k)-q- e* (k) fc") . (51) 

This, in turn, gives [from Eq.(41)] 

F 3 (g 2 ) = 2 [-g + - q 2 h - (M 2 B - q 2 ) h 2 ] (52) 
Similarly, from Eq.(46), we get the relation 

(7 (k, e) \ma a ^b\ B(p)) = -ie aPpc e* p q^f2 [g + - q 2 h x - M 2 B h 2 
Comparison of this equation with Eq.(40) gives 

Fi (V) = (q 2 ) - q 2 h (q 2 ) - M 2 B h 2 (g 2 )] (53) 
Thus, finally we obtain 



Fv (g 2 ) 
F A (?) 



m b + m q 
2 

m b - m q 



{ 9+ {q 2 )-q 2 h 1 (q 2 )-M 2 B h 2 (q 2 )}, (54) 
{g + {q 2 ) ' q'h (q 2 ) - (M 2 - q 2 ) h 2 (g 2 )} . (55) 
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Therefore, the normalization of Fy and Fa at q 2 = is determined by a 
universal form factor (g + (0) — M|/i 2 )- Now the form factor hi does not get 
any contribution from quark triangle graph nor from the pole and therefore 
we shall put it equal to zero. On the other hand, only g + (q 2 ) gets contribution 
from quark A-graph [24], 
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We expect the Ward Identities to hold at low q 2 below the resonance regime 
and as such we use the results obtained from them at q 2 = 0. Thus from 
Eqs. (54 and 55), we obtain 

(m b + m q ) F v (0) = 2g + (0) = (m b - m q ) F A (0) . (57) 

Further, using Eq. (28) in the above Eq. (57) and neglecting terms of the 
order of (A =F m q ^j /M B , we obtain another constraint using Eqs. (38, 39) at 

Now if we restrict to one radial excitation (M 2 = Mi) we obtain from Eq. 

(58) 

(M 2 /M 2 -1)M { ' 

F ( i\ 9+ (°) (ac\\ 

V g > M{l-q 2 /M 2 ) (1 - q 2 /M 2 ) 1 } 

Restoring the subscripts and using the definitions (26) 

2<7+(0) M 2 , 



9b*b~i 



Mb f B «(M 2 B ,jM 2 B ,-l) 
2<?+(0) 



f B (MyM 2 * - i) 



(61) 



while 



jB* A B~y ~ -T7 7 >T (.02) 

Mb fBX (m^/m^-i) 
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Table 1: £?-mesons masses in GeV [30] 







M 


Mi/M 


M 2 /M 


M B 


0" 


5.28 


1.14 


1.24 


M b * 


1- 


5.33 


1.14 


1.24 


M b * 

A 


1+ 


5.71 


1.12 


1.22 



Use g + (0) given in Eq. (56) with Q u = 2/3, namely 



**«>) = f& (63, 



we have the prediction 



9b* s 7 = it t- — -y (64) 

3A [My M 2 B , - l) 



Further 



*V(* 2 ) = ^ ^ V (65) 

^ 2 ) = ^ V (66) 

MB (l-qVM B1 ) (l-qyMy 

This is the final expression for the form factors of our process B — > 7/z//, if 
we restrict to the one radial excitation. We also observe the approximate 
equality Fy(q 2 ) = Fjy(q 2 ) of the form factors which also occur in some other 
models [12, 13]. For numerical work, we shall use 5-meson masses given in 
Table 1 and f B = 0.180 GeV. 

This gives the prediction from Eq.(64) 

9b* Bi = X = 5 - QGeV ~^ (67) 

for A = 5.28 - 4.8 = 0.4 GeV^ 1 [see Eq. (28) and Table 1]. Also, we obtain 
from Eq. (63) 

0+(O) = ^ = 0.15. (68) 
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Further from Eq.(62) 



fB* A B~f 



f B M B * A 2.6 



f B 



A 



6.5 



A 

f B M B * A 



GeV- 1 



(69) 



We now study the effect of the second radial excitation. We go back to 
Eq. (38) and use the constraint (58) to obtain 



*■(«») = 



R^_lUMi_1^ 2 ? 2 i 2g+(0) (-, _ 2 ( 1 , 1 M\ \\ 



(f - g 2 /Mf) (1 - gVMf) (1 - g 2 /M 2 ) 



If we parametrize R as 



2g + (0)1-(1-M 2 /M 2 )A 



M (Mf/M 2 - 1) ' 
where A is a parameter which in principle can be obtained when g B * Bl and 



fB^B'y become known. Then 

2 9+ (0) 1-^(1 + (1 



M' 



M 



q 2 /M 2 ) (1 - q 2 /M 2 ) (1 - q 2 /M 2 ) 



(70) 



For Mi = M 2 the above Eq. (70) reduces to Eq. (60). So the couplings of B 
with 5*7 and B\^i become 

2g+i0)Ml [l " (l - M 2 ,/M 2 ,) (l - M^JMI) A 



9b*b~i 



M B f B * (M 2 ,/M 2 , -l) 



(l - M|./M|.) (l - M 2 B{ /M 2 B .) A] h.QGeV- 1 



1 - ( 1 - M 2 ,/M 2 , ) ( 1 - M B * /M B * ) A 



(71) 
S.hGeV- 1 

(72) 



and the corresponding form factors become 
Fy{q 2 ) 



29 ■+(()) 1 ^ 



1+1-^ 



Mi 



M. 



b (l-g 2 /Ml 5 )(l-g 2 /M 2 ,)(l-g 2 /M 2 ,) 



(73) 
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, ,2\ _ ^9+ (0) _J| \ - <; \ - Jo / / 

(i - q 2 /m 2 a J (l - q 2 /M^ A J (i - q yM 2 A ) 

For numerical values we shall use A = [i.e., Mi = M 2 \ and A = 3 and 
A = 4.8. The second value of A (= 3) corresponds to estimate of g B *By from 
vector meson dominance 

2 f _ 

0b*b 7 = o9b*b p -- E j = 2.76GeV~ 1 
3 

where = v / 2(H) GeV" 1 obtained in [31] and f p -/m p = 205 MeV. 

The third value of A (= 4.8) gives more or less the width for B* — > B*y 
obtained from MI transition in non relativistic quark model (NRQM). These 
values give decay width for B* — > B^y transition 23 keV, 5.5 keV and 0.8 keV 
respectively while MI transition in NRQM predicts it to be 0.9 keV. These 
predictions are testable when above decay width is experimentally measured. 



6 Decay distribution 

The Dalitz plot density 

d 2 T _ d 2 T IB d 2 T sd d 2 T IN T 
dxdy dxdy dxdy dxdy 
= PiB(x,y) + p SD (x,y) + p INT (x,y) (75) 

is a Lorentz invariant which contains the form factors Fy and Fa in the 
following form [19, 20, 22] 

Pib(x,d) = A IB f IB (x,y) 

p SD (x,y) = A SD M 2 B [(F v + F A ) 2 f SD+ (x,y) + (F v -F A ) 2 f SD -(x,y) 
PiNr(x,y) = A INT M B [(F v + F A )f INT +(x,y) + (F v -F A )f INT -(x,y)] 

where 

fiB{x,y) = ( 2 , \ y + Tl r) 

\x z {x + y — 1 — ri) J 
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fsD+(x,y) 
fsD-(x,y) 



(x + y - 1 - ri) ((x + y - 1) (1 - x) - r t ) 
(l-y + r l )((l-x)(l-y)+r l ) 



fiNT- (x, y) 



fiNT+ (x, y) 



{ , \ V \ Tl \ {{l-x){l-x-y)+r l ) 
\x(x + y- l-ri)) 

( < \ — \) ( x2 -(l-x)(l-x-y)- n) 
\x(x + y - 1 - n)) v > 



and 



A 



IB 




The SD + term reaches its maximum at £ = 2/3, y = 1, which corresponds 
to 6ij = 7r. The SD~ term reaches its maximum at x = 2/3, y = 1/3, 
corresponding to 6^ = 0. Indeed, for lepton of maximal energy (y = 1), 
only "right-handed" photons contribute. In this situation, the photon and 
the neutrino must be emitted in the direction opposite to that of the lepton. 
Angular momentum conservation forces the photon spin to be opposite to 
the total lepton spin and the photon helicity has the same sign as that of 
the lepton. Then the photon and the neutrino are emitted parallel. This 
configuration corresponds to a neutrino of maximal energy (E u = E™ ax when 
x + y — 1). In this case, only the "left-handed" photon contributes. When 
x + y = 1, the IB contribution becomes very large: this corresponds to 6*/ 7 = 
0. Consequently, it is very difficult to distinguish experimentally between the 
IB and the SD~ contribution. To summarize, an experiment performed in 
the region 9i y ~ n is essentially sensitive to (F v + Fa) 2 - 

The form factors calculated in Eq. (60) can be expressed in terms of 
dimensionless variable x, 



F(x) 



F(0) 



(76) 



x 1 - (1 - x) I (Mi/M) 



21 ' 
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where x is defined in Eq. (12) and q 2 in Eq. (15). After restoring subscripts, 
the form factors Fy (q 2 ) [Eq. (65)] and Fa (q 2 )[Eq. (66)] can be written as 

Mo) 



F v (x) 



F A (x) 



x 



{l-x)/(M Bl /M B *y 
MO) 



x 



i-{i-x)i(m b1 jm b * 



(77) 
(78) 



where 



'v,a 



(0) 



2fl+ (0) 

Mr 



We use these in Eq. (75) and integrate over x and y in the limit as mentioned 
in Eqs. (16, 17). IB contribution diverges for the minimum value of x, we 
take an arbitrary lower limit for x i.e. x m i n ~ r\ for which the divergence 
problem is cured and the IB part gives some definite value O(10~ 20 ). But as 
the energy of the photon is increased, it approaches zero at x max . Therefore 
in the total decay width, this does not contribute much. The SD part is 
the most dominant part of the decay width which provides almost the whole 
contribution. This part increases initially with increasing x, reaches its peak 
value and then starts decreasing. The INT part of the decay width is an 
increasingly vanishing contribution and can be neglected in comparison to the 
SD part, because it is suppressed by O(10~ 21 ) and becomes flat (approaches 
zero) as x (the photon energy) approaches 1 (its maxima). Therefore, this 
does not contribute fairly to the total decay width of the process. 

In the Fig. 2, differential decay width of the process is plotted against x 
and we see that for our calculations, the peak is shifted to lower value of x as 
compared to those for Eilam et al., [10], Korchemsky et al. [12] and Chelkov 
et al. [13]. So, for the process B — > 7/^ the branching ratios obtained is 

B(B -> = 0.5 x 10- 6 (I = fi) (79) 

This value is for the form factors given in Eqs. (77, 78) which are obtained 
by restricting to the first radial excitation only. Now if we consider the effect 
of second radial excitation the expression for the form factors are given in 
Eqs. (73, 74). The branching ratio thus obtained are 

B(B -> = 0.38 x 10- 6 (/ = /1, A = 3.0) 

B(B -> = 0.32 x 10~ 6 (1 = /1, A = 4.8) 
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for two representative cases of A = 3 and A = 4.8 respectively. These are not 
sensitive to the values of A in contrast to the decay width of B* — > Bj. The 
CLEO collaboration indicate an upper limit on the branching ratio B(B + — > 
7z//e + ) of 2.0 x 10~ 4 at the 90% confidence level [1]. The predicted values are 
within the upper limit provided by CLEO collaboration but differ from those 
predicted in [12, 13], namely (2 — 5) x 10~ 6 and 0.9 x 10~ 6 , respectively. The 
Monte-Carlo simulation results are given in [32] where the upper limit on the 
branching ratio for this process is predicted to be 5.2 x 10~ 5 . 



7 Conclusions 

Preliminary data from the CLEO Collaboration indicate an upper limit on 
the branching ratio B(B + — > 7f/e + ) of 2.0 x 10~ 4 at the 90% confidence level 
[1]. With the better statistics expected from the upcoming B factories, the 
observation and experimental study of this decay could become soon feasible. 
It is therefore of some interest to have a good theoretical control over the 
theoretical uncertainties affecting the relevant matrix elements. 

We have studied B — > 7//// decay using dispersion relations, asymptotic 
behavior of form factors and Ward Identities. The dispersion relation in- 
volves ground state B* and B* A resonances and their radial excitations which 
model contributions from higher states and continuum contribution, which 
is calculated from quark triangle graph. The asymptotic behavior of form 
factors and Ward Identities fix the normalization of the form factors in terms 
of universal function g + (0) at q 2 = and put constraints on the residues. 
Thus in our approach, a parameterization of q 2 dependence of form factors 
is not approximated by single pole contributions. This parameterization 
is dictated by considerations mentioned above and also the coupling con- 
stants of 1~(B*) and 1 + (B\) resonances with photon are predicted if we 
restrict to one radial excitation. By using A = 0.4 GeV -1 we have calculated 
g + (0) = 0.15 and predicted the value of #_b._b 7 = 5.6 GeV -1 (cf. Eq. (67)) 
and fs* A By = 6-5fBM B * A / /b* GeV _1 (cf. Eq. (69)). Taking into account one 
radial excitation the form factors are summarized in Eq. (65, 66). Branch- 
ing ratio for the process is then calculated to be B(B — > 7^/) = 0.5 x 10~ 6 
which lies within the upper limit predicted by CLEO Collaboration at 90% 
confidence level [1]. Then we study the effect of second radial excitation in 
terms of a single parameter A, which in principle is determined once <?_b*_b 7 
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and is* b 7 are known (cf. Eq. (71, 72)). The resulting form factors are 
given in Eqs. (73, 74). By using these form factors the branching ratio is 
B(B -> = 0.38 x 10~ 6 and B(B -> = 0.32 x 10~ 6 for two represen- 
tative cases A = 3.0 and A = 4.8 respectively. These branching ratios are not 
sensitive to the value of A in contrast to radiative coupling constants which 
give respectively B* -> £7 width as 23 keV (A = 0), 5 keV (A = 3.0) and 0.8 
keV (A = 4.8). One can also predict radiative widths of radial excitation in 
terms of B* and B* A radiative widths by using relations (36), (39) and (58). 
The differential decay width versus photon energy is plotted in Fig. 2 to 
compare our results with the existing calculations in the light-cone QCD ap- 
proach [10, 12] and in the instantaneous Bethe-Salpeter approach [13]. The 
results for B — > 71/jZ have been reproduced by using Suddakov resummation 
[12] and have also been shown graphically. In our calculations as well as in 
[10], the position of the peak of differential decay width is shifted towards 
the lower value of photon energy spectrum. This is due to the double pole in 
the form factors. The over all effect of radial excitations is to soften the q 2 - 
behavior of differential decay distribution while in [12] it is due to Suddakov 
resummation. 

Our main inputs have been dispersion relations, asymptotic behavior and 
Ward Identities, all of which have strong theoretical basis and in these aspects 
it differs from other approaches. Our approach is closer to the one followed 
in [18] for B — > -nlv\. Only external parameters involved are fs, resonance 
masses (which are determined in potential models) and (?b*b 7 and /b^b 7 
which are either predicted or on which we have some theoretical information. 
The radiative widths of radial excitations are predicted in terms of the above 
coupling constants. Thus our approach has predictive power and can be 
tested by future experiments. 

The experiments at the B-factories, BaBar at SLAC and Belle at KEK 
(Japan) and the planned hadronic accelerators are capable to measure the 
branching ratio as low as 10~ 8 [33]. 
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Figure Captions: 

1. B — > luff radiative leptonic decay diagrams. 
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2. Differential decay rate versus photon energy x is plotted and com- 
parison is given with various approaches. The solid line (for A = 0), 
dashed-trippledotted line (for A = 3.0) and dotted line (for A = 4.8) are 
our calculation, dash-dot-dot line [10] , dashed line [12] and dash-dotted 
line [13]. The thin-solid line is the Sudekov resummation calculation 
result from Ref. [12]. 
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